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Abstract. The Volterra integral equations of the first kind with piecewise smooth kernel are 
considered. Such equations appear in the theory of optimal control of the evolving systems. 
The existence theorems are proved. The method for constructing approximations of paramet- 
ric families of solutions of such equations is suggested. The parametric family of solutions is 
constructed in terms of a logarithmic power asymptotics. 
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Introduction 

The theory of integral models of evolving systems was first initi- 
ated by V. Glushkov in the 70s of XX century. Readers may refer to 
the bibliorgaphy in HHY961 [ApaOlj IDL951 iMSTTTl IHY031 . Such the- 



ory employs the Volterra integral equations of the first kind where limit 
of integration is time function. The theory and numerical methods of 
such nonclassic equations were studied in the monograph [ |Apa03[ and 
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applied for estimation and optimization of lifetime of electric power 
systems (EPS) components HMST11H . 

In this paper we address the following integral equation: 

J K(t, s)x(s)ds = f(t), < t < T, (1) 
o 

where kernel is defined as follows: 

r ifi(M), o< s < ai (t), 

K u s \ = J K 2(t,s), < s < a 2 (t), (2) 

[ K n (t, s), a n -i(t) < s <t, 
< «i(t) < a 2 {t) < ■ ■ < a n -i(t) < t, \a[(t)\ < 1. 



Functions Ki(t,s),i = l,n,f(t) are continuous and sufficiently 
smooth, /(0) = 0. 

The Glushkov integral model of evolving systems HHY96[|MST11[ 
IHY03D is the special case of the Volterra integral equation (OQ) - © 
when all the kernels except Ki(t, s) are zeros. 

First results in studies of the Volterra equations with discontinuous 
kernels were formulated by G.C. Evans HEvalOl in the beginning of XX 
century. Results in the spectral theory of integral operators with discon- 
tinuous kernels were obtained by A. P. Khromov in HKhr06l . Asymp- 
totic approximations of solutions of the Volterra equations of the first 
kind with analytical kernel K(t, s) were constracted by NA. Magnit- 
sky HMag83| l. 

It is to be noted that solutions of the equations (OQ) can have an arbi- 
trary constants and can be unbounded for t — > 0. For example, if 

J ^' s) = \ -l,f/2< a <f, (3) 

f(t) = t, then equation (OQ) has the solution x(t) — c — p|, where c is 
constant. 



The Volterra Models 3 



In this paper we employ results of the papers [Mag83[rSS06[|ST09l 



ISST07H in order to formulate the algorithm for construction of the con- 
tinuous solutions of equation (Q~|) for < t < T in the following form: 



N 



x (t) = ^XiQnt)? + t N u(t). (4) 



i=0 



Coefficients Xj(lnt) are constructed as polynomials on powers of hit 
and they may depend on certain number of arbitrary constants. N de- 
fines the necessary smoothness of the functions Ki(t, s), /(£). In this 
paper we propose an algorithm for construction of the function u(t) in 
representation of the desired solution © based on successive approxi- 
mations method which is uniformly converge on [0, T]. It is to be noted, 
that logarithmic-power asymptotics have been efficiently employed for 
solution of integral and differential equations in irregular cases HST09M 
-ESID, Pig83] |. 



The paper is organized as follows. In Section 1 after the problem 
statement we introduce the structure of solution and prove the existence 
theorem. The method for the asymptotic approximations construction 
is suggested in Section 2. The main theorem is formulated and proved 
also in Section 2. Finally, concluding remarks are given. 



1. The structure of solutions and existence theorem 

For sake of clarity let us suppose that a^t) = a^, where = a < 
ax < ■ ■ ■ < a n — 1. Let us introduce the condition 
A. K n (t, t) ^ for t E [0, T] and N is selected to fulfill the following 
inequality 

n , 

max \K n (t, t)]- 1 ( <*l +N \Ki{t, (5) 
i=i ^ 

+^+ JV |# i (t,ai_ 1 t)M <l + q, 

where q < 1, ao = 0, a n = 1. Condition A is fulfilled for large enough 
N since on G (0, 1) for i = 1, n — 1 
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Lemma 1.1 Let condition A be fulfilled, let all the functions 
Ki(t,s),i = l,n be differentiable wrt t and continuous wrt s. Then 
the homogenius equation 

t 

J K(t,s)s N u(s)ds = (6) 
o 

has the trivial solution in the space C[o,t]- 



Proof 1.1 Let us differentiate the equation © and take into account 
(|2|). Then we get an equivalent integral-functional equation 



K n (t,t) \t 

% <Xi-\t 



where 



n , 

Lu = ( a] +N Ki(t, ait)u(aiit)- 



-al+?Ki(t, a i _ 1 t)u(a i .. 1 t) {K n {t, t))- 1 



Due to the condition A in the space C^ti we have the following esti- 
mate 

\\Lu — u\\ < 

Therefore, according to the theorem on the inverse operator ( l\Tre07\I . 
p. 134), and because of the inequalities = ao < ai < ■ ■ ■ < a n = 1 
there exists the following bounded inverse operator L^ 1 £ £(C[ 0j t] — > 

C[0,T]) 

II^II^T^- (8) 
1 - q 

and equation © can be reduced as follows 

at 

u{t) = -L~^ I ^^(s/t) N u(s)ds = Au, (9) 
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where < t < T. Let us introduce the following equivalent norm 
I lull = max e~ lt \u(t)\, I > in the space Qoti- In this norm the 

0<t<T 1,1 

inequality (IS]) remains correct and for sufficiently large I an operator 
A will be contracting since \\A\\ < y^m(Z), where m{l) — > for 
I — > +oo. Therefore homogeneous equation (|9]) has trivial solution. 

Corollary 1.1 Let all the condition of Lemma 1 be fulfilled, g(t) G 

CjgL , \g'(f)\ = o(t N ) for t — > +0. Then inhomogeneous equation 

t 

J K(t, s)s N u(s)ds = g(t) has the unique solution, and u(t) — )■ for 
o 

t ->■ +0. 

Proof is trivial since differentiation of this equation leads to the 
equivalent equation 

u(t) = Au + r N L~ l g'{t) (10) 
with contracting operator A and continuous free function. 

Theorem 1.1 Let in the space of continuous on (0, T] functions which 
have the finite limit fort — > +0 (briefly, in class C(o,t]J exists the func- 
tion x N (t) such as fort —> +0 

t f 

(- J K(t,s)x N (s)ds + f(t^j =o(t N ), 



f(t) G C[q T ], /(0) =0. Then equation (E]) has the following solution 

x(t) =x N (t)+t N u(t) (11) 

in class C( ,t] • Here function u(t) G C[ 0i t] , — > 0/or t — )■ +0 and it 
can be uniquely constructed with successive approximations method. 

Proof 1.2 Proof follows from the corrolary 1. Indeed, with 4771) we can 
rewrite the equation (UJ) as follows 

t 

J Kit, s)s N u(s)ds = g(t), (12) 
o 
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where function g(t) is following 

t 

g(t) = - J K(t,s)x N (s)ds + f(t) (13) 
o 

and it satisfies the condition of the corollary 1. Therefore in ([77]) the 
function u(t) can be uniquely constructed with successive approxima- 
tions from the equation rfTOl) using arbitrary initial condition. 

Definition 1.1 The equation ([72]) with right hand side ([72]) we call reg- 
ularization of the equation (UJ), and function x N (t) as asymptotic ap- 
proximation of solution ( 1771) of the equation ([/J. 



It is to be noted that one could numerically find the function u(t) by 
solution of the equation (IT2T) based on well known numerical quadra- 
ture schemes (see e.g. the bibliography in the monograph [Apa03|). 
The method of constructing asymptotic approximations x N (t) in the 
solution (fT2l) we will study below in the Section 2. 



2. The method of asymptotic approximations construction 

Let us suppose that along with the condition A the condition B be 
fulfilled. Functions Ki(t,s),i = l,n, f(t) are N + 1 times differen- 
tiable in the neighborhood of zero, where N is selected according to the 
condition A. We introduce an auxiliary algebraic equation wrt j E N 

n 

± £ K t (0, 0)(a] + > - a]+[) = (14) 
i=i 

and name it as characteristic equation of the integral equation (Q])- Since 
f(0) = 0, then equation 

n 

'^(a i K i (t,a i t)x(a i t)) - a^K^t, a i _ 1 t)x(a i „ 1 t))) + 

i=l 

n I 

+ K' i (t,s)x(s)ds = f(t) 

i=l ^ , 

ati-it 
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is equivalent to the equation (Q])- We will look for the asymptoti- 
cal approximation of it's solution as following polynomial x N (t) = 

N 

Xj(\nt)P . Based on the method of undetermined coefficients, and 

3=0 

taking into account the inequalities = ao < cm < • • ■ < a n = 1 we 
construct the recursive sequence of difference equations wrt the coeffi- 
cients xj(z)(z = Int) as follows: 

K n {0, 0)x 3 -,(z) + (#*(0> 0)- 

i=l 

- K i+1 (0, 0))xj(z + a,i) = Mj(x , Xj-i), (15) 

where j = 0,N, a,i = In c^, i = 1, n — 1, M = /'(0). 
Here we follow ( HGel59L p. 330) and we seek the solution of the ho- 
mogeneous difference equations in the form of x = \ z . Substitution of 
the function \ z into the homogenius difference equations leads to iV + 1 
equations for difference equations (PT51) : 

n-1 

V S {X) = K n (0,0) +^^'(^(0,0)- (16) 

i=l 

K l+1 (0,0))\ a > = 0, j = 0~N. 

Therefore we have 

Property 1. jth equation H6\) has the root A = 1 if and only if j satis- 
fies characteristic equation rf7?l) of the integral equation (0). Moreover, 
multiplicity of the root j of the equation rf7?l) is equal to rj iff 

n 

L{3) = ^^(0,0)K 1+ ^ - a£}) = 0, (17) 

i=l 

n-1 

«i +i (^(0, 0) - K i+1 (0, 0))a\ = 0,1 = 1,^-1, (18) 

n-1 

5>J +3 '(tf<(0,0) - K i+1 (0,0))a? jt 0, 

;=i 
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where a = 0,a n = l,a = 0,a n = 0, a« = lnajj, i = 1, n — 1, and 
multiplicity r, < n — 1. 

Proof follows from the equality 



and from the structure of the equations (fl4l) . (PT6l) . 

If we suppose that for certain j multiplicity r\, > n, then K\(0, 0) = 
if 2 (0,0) = ... ••■ = AT n _i(0,0) = K n (0,0) due to CI, since 
det |K||u=t^ 7^ 0. But due to CLU) -K" n (0,0) = 0, which contradicts 
A. 

Under the conditions A, B there are two cases. 

2.1. Regular case 

Let L(j) ^ 0, j G N. Then A = 1 does not satisfy any of the equa- 
tions in the sequence (PT6l) . All the coefficients Xj of the asymptotics 



x N = Xit % can be determined uniquely with method of undetermined 



coefficients and do not depend upon In t. 
Therefore we have the following theorem 

Theorem 2.1 Let the conditions A, B and ^ 0,j 6 N be fulfilled. 

N 

Then equation (UJ has in C^ti the solution x(t) = + t N u(t), 

i=0 

where Xi are determined uniquely with method of undetermined coef- 
ficients and function u(t) is uniquely constructed (numerically or with 
successive approximations from equation 4721)). 

2.2. Irregular case 

Let L(j) = only for j e {ji, . . . , j k } C {0, 1, . . . , N} and multi- 
plicity of the root A = 1 for the corresponding characteristic equation is 
r,j. Let in jth difference equation (PT5T) right hand side Mj(z) appears to 



n-l 



n 




(19) 



N 
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be polynomial from z of the order rij > 0. Then in irregular case, i.e. for 
Tj > 1, based on ( HGel59L p. 338) particular solution of the jth equa- 

rij+rj 

tion (fT51) we have to search as following polynomial x(z) = Yl c i z% ■ 

i=rj 

Coefficients q of this polynomial can be sequentially calculated by the 
method of undetermined coefficients starting from c nj+r .. Coefficient 
Xj (z) of desired asymptotical approximation x N in this case is as fol- 
lows 

Xj(z) = c + c\z H V c rj _ 1 + x(z). 

In irregular case when rj > 1, constants c , . . . , c Tj _\ remain arbitrary 
since functions z\ i = 0, 1, . . . , rj— 1 satisfy jth homogenius difference 
equation corresponding to (IT3T) . 

In applications one could use the Propetry 1, and coefficient Xj(z) 

n j+ r j 

in irregular case directly as polynomial Yl °i z \ where c nj+rj , c 

i=0 

are determined sequentially using method of undetermined coefficients. 
And c rj _i, . . . , Co remains arbitrary. Therefore in irregular case when 
L(j) =0 for some j new arbitrary constants rj appear in determination 
of the coefficient Xj(z). Order of the polynomial Xj(z) on the value of 
the multiplicity of rj or root A = 1 of the jth equation (PT6l) becomes 
greater than order rij of the right hand side of the corresponding equa- 
tion (fT3T) . i.e. of the order of the polynomial Mj(z). 

Therefore we have the following theorem: 

Theorem 2.2 Let the conditions A, B are fulfilled. Let characteristic 
equation L(j) = of integral equation ([/]) has exactly k natural roots 
{jii ■ ■ ■ ,jk}- And let the root A = 1 of jth equation rf73l) has multiplicity 
Tj. Then equation ((7]) has the following solution in C( 0j t] 

N 

x = ^ Xi (lnt)f + t N u(t), (20) 

i=0 

which depends on p = r± + • • ■ + arbitrary constants. Moreover, 
coefficints Xi of the asymptotic approximation x N (t) are polynomials 
from Int. Function u{t) can be constructed by successive approxima- 
tions which converge uniformly for t G [0, T], or numerically from rf72l) . 
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Remark 2.1 If L(0) = 0, then in solution rf20l) xq = const + a hit, 
where a is the defined constant. Therefore in this case x(t) G 
C(o,t], ^lim^t) = oo. 

Remark 2.2 These results can be generalized if in the equation §Tj) for 

a*-i(t) < s < K(t, s) = s) w/zere a^O) = 0, < a' (0) < 
ai(0) < ••■ < <(0) < l,z = l,..., n. 

Conclusion 

Our method allow studies of the Volterra integral equations of the 
first kind when kernals are discontinuous operator-functions operating 
in Banach spaces. Our work naturally complements the theory of inte- 
gral models of evolving systems. 
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